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ABSTRACT 

We derive a variational principle for the dynamical stability of a cluster as a 
gas sphere in a box. Newtonian clusters are always dynamically stable and, for 
relativistic clusters, the relation between dynamical and thermodynamical insta- 
bilities is analyzed. The boundaries between dynamically and thermodynamically 
stable and unstable models are found numerically for relativistic stellar systems 
with different cut off parameters. A criterion based on binding energy curve is 
used for determination of the boundary of dynamical stability. 
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1. Introduction 

It is well known that Newtonian stellar clusters with effective adiabatic power 7 = 5/3 
are always dynamically stable (see Zel'dovich & Novikov 1971). An isolated cluster, in which 
temperature tends to a constant value all over the radius, suffers gravothermal catastrophe 
(see Antonov 1962; Lynden-Bell & Wood 1968) at which any finite object evolves into a 
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model with highly concentrated core and very extended envelope, where radius and central 
density tend to infinity. If we remove a demand of constant temperature and consider a rapid 
dynamical (adiabatic) perturbation, then the cluster does not react drastically and returns 
always to its mechanical equilibrium. 

An interesting analysis of dynamical and thermodynamical instability of stellar clusters 
made by Chavanis (2002a, 2002b) contains some sections that show this problem is still not 
quite clear, instead of several publications devoted to this topic (see Merafina 1999; Lightman 

6 Shapiro 1978). 

In section 2 we derive in Newtonian gravity the variational principle for investigation 
of dynamical stability of a stellar cluster or a gas in a spherical box and discuss properties 
of trial functions which may be used for a stability analysis. For Newtonian gravity, we find 
conditions 7 > 4/3 for stability of an extended cluster with Pe/Po ^ 1 ^-nd the condition 

7 > for stability of a very hot body in a box at P/p ^ GM / R with almost constant 
pressure P and density p. The pressure along the adiabate is supposed to follow the relation 

P = fcpT. 

Relations between dynamical and thermodynamical stability of relativistic clusters are 
analyzed in sections 3 and 4. Numerical results about stability analysis of relativistic stellar 
clusters with different cut off parameters are represented in section 5. 

The oscillatory behavior of the mass M of the cluster as a function of a central den- 
sity poi at fixed temperature T, indicates increasing number of thermodynamically unstable 
modes. Similar oscillatory dependence of M(po), at fixed parameter Wq, shows (approxi- 
mately) increasing number of dynamically unstable modes. Dynamically stable models with 
arbitrarily large central redshifts Zc exist only at temperature T < 0.06. All such models are 
thermodynamically unstable. 



2. Newtonian clusters and gas spheres 

Let us consider development of dynamical perturbations in the cluster, where character- 
istic time is usually much shorter than any other time, including the time of energy exchange 
with an outer thermostat. In dynamical (rapid) perturbations, where local entropy is con- 
served and there is no time to smooth the temperature, the relations between 5T, 5p and 
SP follow adiabatic relations 

SP _5 5p _5ST 

~¥~Z~p~ 2 7^ ' ^ ^ 
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so the pressure perturbation, where in mechanical equihbrium we have 

dP kT dp 



should be taken as 



dr rUs dr 



(2) 



5P ^ ^-Sp 
P 



(where 7 = 5/3) . 



(3) 



The dynamical stability analysis (Chandrasekhar 1964; Bisnovatyi-Kogan 2001) is usually 
performed for a star with zero density and pressure on the boundary. For a cluster in 
a box with nonzero p and P at the edge, the boundary conditions are different, but the 
expression for frequency, following from the variational principle, remains the same. It is 
easy to show, using equations in Lagrangian coordinates in which dynamical equations at 
spherical symmetry are written as 



dv 1 dP Gm 
dt p dr 



dr 
di 



(4) 



(5) 



dr 



dm Anpr"^ ' 

that, for linear perturbations of the static model with Sr ~ e**^*, we obtain, using Eq. (3), 
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The equation for a perturbation 5r{m), taking into account Eqs. (7), (8) and the equilibrium 
equations, is written as 



AdP5r 'yd 

-a dr H ; 

p dr r pdr 



r dr 



= 



(10) 



The Eq. (10) for a cluster in a box should be solved at boundary conditions 



Sr{0) = 6r{M) = , r(M) = R , 



with a variable outer density and pressure 6p{M) ^ 0, 5P{M) ^ 0. This is different from 
the boundary conditions of stellar oscillations which, in adiabatic approximation, are 



5r(0) = 



5p{M) 



^5r d5r\ 
^ r dr J 



(12) 



0, 



m=M 



with a variable outer boundary radius 6r{M) ^ 0. 



Instead of solving Eq. (10), let us derive a variational principle, which gives a possibility 
of a simple stability analysis. Integrating Eq. (10) over the mass of the cluster, after mul- 
tiplying by 6r, and accepting the normalization of the linear perturbation function in the 
form 



M 



5r^dm — A , 



(13) 



we obtain from Eq. (10), after partial integration with boundary conditions (11), the follow- 
ing expression for the squared frequency 
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This is exactly the same expression which takes place for a star with another boundary 
conditions (12). 
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Variational principle over-estimates the values of for different trial functions 5r(m), 
giving the minimal exact value for the eigenfunction of the oscillations. So, variational 
principle may prove the existence of instability, but only approximately permits to make a 
judgement about the stability of the system. It is important to use only those trial functions 
which satisfy the boundary conditions of the eigefunction. 

It is known for stars with boundary conditions (12), that the linear trial function gives 
almost an exact result for stability boundary. We have from Eq. (14) for — ar 



that corresponds to stability boundary at 7 = 4/3. This boundary is the exact value for stars 
with constant 7 because at 7 = 4/3 the trial function 5r — ar is also an exact eigenfunction 
(Zel'dovich & Novikov 1971). 

For a cluster in the box the linear eigenfunction is not vahd, because it does not sat- 
isfy outer boundary conditions. Nevertheless, for clusters with low ratio Pe/po (where pe 
and po are the external and central density, respectively) the exact fulfilment of the outer 
boundary condition is unimportant and, at Pe/po — 0, the condition (15) is approximately 
valid also as a criterion for dynamical stability of a cluster in a box. On the other hand, the 
dynamic response of a stellar cluster with quasi- maxwcllian distribution function to global 
radial perturbations is similar to the response for adiabatic perturbations in a star with 
the same adiabatic index, which is 7 = 5/3 in the nonrelativistic cluster. In fact, for non- 
relativistic spherical stellar clusters, we obtain dynamical stability because also barotropic 
stars with the same density distribution function are dynamically stable (Antonov 1960). 
This result leads in many cases to choose a parallel treatment in stability analysis for these 
physically different systems (Binney & Tremaine 1987). Similar correspondence exists even 
in the investigation of thermodynamical stabihty (Bettwieser & Sugimoto 1985). Therefore, 
nonrelativistic clusters with 7 = 5/3 in the process of gravithcrmal catastrophe arc, always, 
dynamically stable. Chavanis (2002a) considered perturbations at constant temperature and 
obtained dynamic instability in this case. We should stress that perturbations developing in 
the dynamical time do not preserve constant temperature, which could be reached only by 
contact with an external thermostat. Therefore the increment of "dynamic" instability at 
constant temperature is determined by the time of heat exchange between the cluster and 
thermostat, but not by the time of dynamical processes inside the clusters. 

For clusters in a smaller box the outer boundary condition plays more important role 
and, for small boxes, where the gravity is less important (in the limit P — constant), using 




(15) 
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the trial function 5r — ar gives an highly inadequate result. Using a trial function with 

a correct boundary condition is necessary here, what, as evidently expected, reduce the 
demand on 7 for stability of the object in the box. Therefore, only a relativistic gas in the 
box may become dynamically unstable, because the Newtonian clusters have always 7 = 5/3. 
Using a trial function for a box with radius R and constant pressure Pq as 



Sr — ar{l — r/ Ft) 



(16) 



we get from Eq. (14) 
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where 7 is constant all over the configuration. 

The stability criterion 7 > corresponds to the dynamical stability of a weakly gravi- 
tating gas in the box under the condition of a high pressure Po/p ^ GM/R. The criterion 
7 > remains valid for a more general trial function 6r = Q;r"(l — r/R), for any value 
of n > 0. We see that all instabihties of the Newtonian isothermal clusters with different 
boundary conditions considered by Antonov (1962), Lynden-Bell & Wood (1968) and Cha- 
vanis (2002a) are related only to thermodynamical instabilities, being the clusters always 
dynamically stable. 



3. Relativistic models of selfgravitating systems 

Gaseous systems in a box 

The relation between dynamical and thermodynamical stability is much tenser in the 
relativistic case, when increased gravitational force may lead to dynamical instability for any 
equation of state. 

The analysis made by Chavanis (2002b) demonstrates the development of dynamical 
instability for stars with ultrarelativistic equation of state in the box at increasing box size. 
We should clearly distinguish between two limits of the ultrarelativistic equation of state 
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qe 



(18) 



considered by Chavanis. 

The first one corresponds to the deep interior of the neutron star where, in conditions 
of ultrarelativistic degeneracy and strong nuclear forces, we have Eq. (18) with 1/13 < q <1 
(see Ambartsumyan & Saakyan 1961; Zel'dovich 1962a). The parameter characterizing the 
dynamical stability is the adiabatic index (see Harrison et al. 1965) 



In Newtonian limit we have always 7cr = 4/3 but, when the effects of general relativity 
make the gravity stronger, 7cr becomes larger than 4/3 (see Kaplan 1949; Chandrasekhar 
1964; Merafina & Ruffini 1989). Then, isolated ultrarelativistic stars for which P ~ p(? 
are always dynamically unstable, while stable existence of such configurations is possible, in 
principle, only inside a box with fixed radius or fixed external pressure. Clearly, in conditions 
of zero temperature there is no sense to discuss the thermodynamical stability of the system. 
Therefore, the analysis made in Section 3 and 4 by Chavanis (2002b) relates only to the 
dynamical stability of such star (see also Yabushita 1974). 

The loss of stability, dynamical as well as thermodynamical, may be found from the 
linear series of equilibrium models, where the extremum of mass in the appropriate curve 
M(po) determines the appearance or disappearance of unstable mode. For dynamical stellar 
stability this "static" criterion was formulated by Zel'dovich (1963) and, for thermodynamical 
instability, was used by Lynden-Bell & Wood (1968). Analysis of Chavanis (2002b) has 
shown that dynamical instability of relativistic sphere in a box with fixed radius happens 
exactly in the first maximum of the function M(po)- Contrary to that, Yabushita (1974), 
who investigated dynamical stability of gas spheres in a box with fixed external pressure in 
general relativity (GR), had found the loss of stability in a point well before the maximum 
of the curve M(po)- It is important to note that not every kind of linear series of models 
may be used for investigation of dynamical stability. In the case of an isolated star with 
a zero boundary pressure it is necessary to have a fixed distribution of specific entropy 




(19) 



where, for dynamically stable configurations, we have 



7 > 7, 



cr 



(20) 
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S{m/M) along the series of models M(po), where the first maximum of M{po) denotes the 
loss of stability relative to the global mode without nodes (sec Zel'dovich 1963). Analysis of 
dynamical stability of stars in a box with fixed radius and fixed external pressure had shown 
that only the first maximum of the curve M(po)|re corresponds to the loss of dynamical 
stability, but maximum on the curve M{po)\p^ is situated after the point of the loss of 
dynamical stability. The spheres in a box with constant external pressure considered by 
Yabushita (1974) are not isolated from the surrounding medium, which produces a work 
when the box is changing its volume. So the comparison of two models with equal masses 
and different central densities has no sense because these models have different energies 
due to different volumes of the box. Therefore, the first maximum on the sequence with 
constant Pi,, considered by Yabushita does not correspond to the onset of instability in spite 
of constant specific entropy of the matter considered along the series. The coincidence of the 
maximum of the curve M(po)|pf, with the point of the loss of dynamical stability happens 
only at g = 0, that formally corresponds to zero external pressure and no external work. 

In order to have a possibility to judge about the onset of dynamical instability from the 
linear series of models in presence of constant external pressure Pb at the stellar boundary 
with a nondimensional radius (usual Emden coordinates), we must take into account the 
work of the external pressure and use the function 



where V is the total volume of the model, instead of M{po). It is easy to show, using 
Yabushita results and notations, that for these models the extrema of the functions S{pq)\p^ 
or Ph{C,b)\£ coincide exactly with the point of the loss of stability for different modes, with 
increasing po or ^i, (see also Chavanis 2003). 



The second limit to which relativistic analysis of Chavanis (2002b) is apphed concerns 
"isothermal" relativistic star clusters, with a local constant temperature (Bisnovatyi-Kogan 
& Zel'dovich 1969; Bisnovatyi-Kogan & Thorne 1970). Oppenheimer-Volkoff equations de- 
scribing the equilibrium of a relativistic cluster are given by 



S{po)^M + 



PbV 



(21) 



Star dusters with cutoff 
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with P(0) = Pq and Mr{0) — 0. Pressure P and total energy density pc^ are expressed as 
integrals in momentum space with the distribution function (29), is the mass inside the 
lagrangian radius r (Bisnovatyi-Kogan et al. 1993, 1998). The Schwarzschild-type metric 
was chosen 

ds^^e'^c^dt^-e^dr^-r^de^ + sin^edip^) , (23) 
where the metric coefficients are defined by the expressions 



e'' = exp ( 2 / „ / dr 



dP/dr 
P + pc^ 



(24) 



In our calculations we have used the variable W instead of P, which leads to equilibrium 
equation (Merafina & Ruffini 1989) 

dW G fl- (3W\ MrC^ + 47iPr^ 

(25) 



dr c^\l3J r{rc^ - 2GMr) ' 

with the condition for integration 1^(0) = Wq. Here (3 — T^lrxK? and W is defined by (27). 
Therefore, each model is uniquely determined by choosing the parameters Wq and T (or /3). 

These clusters are not in exact thermo dynamical equilibrium, so the investigation of 
the behavior of the linear series of models of such clusters in a box does not give, rigorously 
speaking, correct results about neither thermodynamical nor dynamical stability. Never- 
theless at T < mc^, when the cluster is almost nonrelativistic and gravitational potential 
<^ c^, the local temperature is almost constant all over the cluster and therefore the anal- 
ysis of Chavanis (2002b) can give a valid presentation about thermodynamical stability of 
such system. In order to investigate dynamical instability of such clusters, the equations 
for small perturbations should be solved or variational principle may be used: this method, 
which is more complicated in the relativistic case than in the Newtonian one, was derived, for 
relativistic stars, by Chandrasekhar (1964) and Harrison et al. (1965); for relativistic clus- 
ters, the variational principle derived by Ipser & Thorne (1968) and completed by Fackerell 
(1970) was used by Bisnovatyi-Kogan & Thorne (1970) for the investigation of dynamical 
stabihty of the "isothermal" clusters. 
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As it was noted by Lynden-Bell & Wood (1968 and references therein), the behavior of 
a star cluster in a box is related, in some sense, to the isolated cluster with cutoff in the 
energy distribution. Relativistic Maxwellian clusters with cutoff have been first introduced 
by Zel'dovich & Podurets (1965) and their dynamical stability was studied by Bisnovatyi- 
Kogan et al. (1993, 1998). Due to relativistic gravity, the loss of dynamical stability happens 
at (7) < 7cr and the structure of the marginally stable model strongly depends on the cutoff 
parameter. In the relativistic cluster, the lower adiabatic index 7 < 5/3 is connected with 
relativistic motion of an ideal nondegenerate gas of stars: it is decreasing with the increase 
of the temperature. Contrarily to that, a fully degenerate highly non ideal nuclear matter is 
present in the interior of a neutron star which may have even larger adiabatic index 7 > 5/3. 

An analogy exists between the cutoff parameter Wq used by Merafina & Ruffini (1989) 
and the quantity Vi — m{(po — (Pe)/T used by Lynden-Bell & Wood (1968), where (po and 
(Pe are the Newtonian gravitational potential in the center and at the edge of the cluster 
respectively. In a full equilibrium cluster in the box with the density distribution p = 
poexp [m(po — m(p{r)] /T, the quantity —vi represents the logarithm of the ratio of densities 
in the center and at the edge of the cluster 



-.,^^^^--^°^^ln^. (26) 

Pe 

Note that density is a finite value at the outer boundary. The ratio po/Pe is also called 
"density contrast" . The value Wq is defined as 



Wo = ^ , (27) 

where = Te~'^^^^^'^ is the local thermodynamical temperature and the constant T is the 
temperature for an infinitely-remote observer. The metric coefficient i/(r) has an usual 
meaning as in the Schwarzschild metric (Landau & Lifshitz 1962). It is easy to show (see 
Bisnovatyi-Kogan et al. 1993) that, in Newtonian limit, Wq formally reduces to —Vi in the 
first relation to the right side of Eq. (26) but, because of deviation from thermodynamical 
equilibrium implied by the cutoff, the outer edge of the cluster corresponds to zero density. 

Thermodynamical stability of the isothermal cluster in Newtonian gravity with a trun- 
cated distribution function may be characterized by the curve Ei,{Wo) at constant T, where 
its maximum denotes the loss of thermodynamical stability, Ei, — —Em/MT is a nondimen- 

sional specific binding energy of the cluster, E = T + U is the total energy of the cluster 
(kinetic and gravitational energy). The value Ere/GM"^, plotted in Lynden-Bell & Wood 
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paper (1968) for an isolated equilibrium cluster in a box as a functions of —vi at constant 

Te, has the same meaning and characterizes the thermodynamical stability of such cluster. 
The plot GM/TfT as a functions of — fi, also at fixed M and rg, going along the sequence of 
models with varying T determines, indeed, the thermodynamical stability of the cluster in 
a thermal bath. The loss of stability in such a cluster (at —Vi = 3.47, with density contrast 
Po/Pe = 32.125) happens before that of the isolated one (at —vi = 6.55, with density contrast 
Po/Pe = 708.61). In the isolated cluster perturbations are developed at constant energy but, 
for the perturbed cluster in the bath, the temperature is preserved, so in the last case the 
cluster occurs to be less stable. 

The definition of the thermodynamical stability or instability of an open cluster with a 
cutoff and zero density at the edge is, generically speaking, senseless, because all such clusters 
are "thermodynamically unstable" . The relaxation in these systems leads to approaching a 
local Maxwellian distribution function without cutoff, equivalent to the formation of a ther- 
modynamically unstable isothermal gas sphere. Existence of maxima on the curve Eh^T{Wo) 
or, equivalently, on the curve Eb^T{po), where henceforth Ej, will simply indicate the specific 
binding energy E^/N, may signify, instead, the appearance of an additional "thermodynam- 
ically" unstable mode which is developed without an increase of cutoff parameter and only 
due to a global fluctuation of the parameter T. In reality, thermodynamical instability of 
both types is governed by the so-called "two-body relaxation time" (Ambartsumyan 1938; 
Spitzer 1940) 



where m is the star mass, R is the radius of the cluster expressed in parsec and is the 
total number of stars in the cluster, and therefore, in this sense, all dynamically stable stellar 
clusters in space are thermodynamically unstable. On the other hand, the local value of 
is proportional to 1/n, where n is the number density of stars, and the influence of cutoff 
is strong near the outer boundary, where local Tb has a maximal value. Moreover, after the 
point of the loss of thermodynamical stabihty, the instabihty begins to develop also in the 
central regions, where local has a minimal value, much lower than at the edge of the 
cluster. Therefore, the loss of thermodynamical stabihty determined by the curve Ei,,t{Wo) 
or £'6,t(po) is important for the cluster with a cutoff almost as well as for the cluster in the 
box. 

The analogy between our truncated clusters and isothermal cluster in the box appears 
if we surround the truncated cluster by the box and let them relax to thermal equilibrium. 
In this sense we exclude the global thermodynamical instability of the open cluster and 




(28) 
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the curve £'6,r(po) or Efj^riWo) gives informations about thermodynamical stability of the 
cluster in a box with a radius equal to the radius of the open truncated cluster. Comparison 
of Newtonian curve of specific binding energy E}JN (see Fig. 1) for open clusters considered 
in the paper of Bisnovatyi-Kogan ct al. (1998) with the corresponding one for clusters 
in a box, indicated in Fig. 2 of the paper of Lynden-Bell & Wood (1968), shows a good 
correspondence between first extrema of these curves which lay at —v\ = 6.55, for clusters in 
a box, and at Wq — 6.42, for open clusters with truncated Maxwellian distribution function. 
This similarity may be seen also from the comparison of these curves as a whole, where 
approximate coincidence of two subsequent extrema is visible. 

4. Relativistic and nonrelativistic oscillations in the curve Eb{po) 
Stability criteria for relativistic models 

The curve Eh^po) characterizing the stability of the star shows a strikingly similar os- 
cillatory behavior for dynamical and thermodynamical types of instability. This oscillatory 
behavior was first analyzed by Dmitriyev & Khohn (1963) in the curve M(po) for cold neu- 
tron stars. They have shown that the star loses its dynamical stability relative to the global 
contraction in correspondence to the maximum of the curve M(po) and each new maximum 
and minimum leads to appearance of a new unstable mode (sec also Misner and Zapolsky 
1964). The detailed analysis of this curve was done by Harrison et al. (1965), where the 
dependence M[R) was also used: this choice is even more useful for analysis of instability, 
because the behavior of the spiral curve M{R) permits to distinguish unambiguously be- 
tween the appearance of a new mode of instability or the removal of the unstable mode, 
both of which may happen in the extremum of the curve M(po)- The analysis of Chavanis 
(2002b), made for a simple relativistic equation of state P — qe in a, box, is very similar 
to the considerations of Dmitriyev & Kholin (1963) and Harrison et al. (1965) made for 
neutron stars. 

The dependence of the total mass M of the star on its total barion number Nj,, which 
has always a positive derivative and shows an angle on the curve M{Nh) corresponding to 
the extremum of the curve M{pq) obtained by Chavanis (2002b) for a star in the box, was 
also demonstrated by Zel'dovich (1962b) for a cold neutron star. 

Similar oscillations for the thermodynamical instability have been obtained by Lynden- 
Bell & Wood for an isothermal cluster in the box. In the paper of Bisnovatyi-Kogan et al. 
(1998) the oscillatory behavior was found in the curve Mr(po) with a related spiral behavior 
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of the curve Mr (a), where the parameter a determines the cutoff of the distribution function^ 

E<mc^ - aT/2 , f ~ e"-^/^ . (29) 

This oscillatory (spiral) behavior, indicated by Bisnovatyi-Kogan et al. (1998) for different 
fixed values of temperature T (Newtonian case corresponds to T — > 0), shows appearance of 
modes of "thermodynamical" instability of clusters with a cutoff, relative to perturbations 
with heat redistribution inside the cluster and change of the cutoff parameter a. More 
detailed curves Mt{po) and Mria) are shown in Figs. 2 and 3, where results are represented 
in nondimensional coordinates from Bisnovatyi-Kogan et al. (1998) with using the same 
calculation scheme. Note that criterion for thermodynamical stability based on the curve 
Mt{po) works sufficiently well only for large T > 0.1. For smaller temperatures the criterion 
based on the curve Eh^T{po) should be used like in Newtonian clusters (see next section). 

Several approximate criteria had been suggested by Bisnovatyi-Kogan et al. (1993) for 
investigation of dynamical stability of rclativistic clusters. All these criteria work almost 
equally well for moderate values of a < 1.5. At larger a, the criterion based on the evalua- 
tion of extremum of the curve Ma{po) is not appropriate because of appearance of increasing 
number of loops (see Fig. 2) which are connected with multiple intersections of the vertical 
line a — constant with the spiral curves of Fig. 3, corresponding to very low (zero) tem- 
perature. The first loop appears at a ~ 1.5, at 1-9 < a < 1.5 there are two loops and, in 
general, there is an even number of loops, except boundaries at a ~ 1.5 , 1.9 , 2.0 , ... , 2.02 
at which there is an odd number of loops Ni = 3 , 5 , ... , oo. At a < ttoo — 2.02 there 
is a curve Maipo) going to infinity and at a > «oo all curves Ma{po) are represented only 
by loops at finite densities. Another important value is — 2.87 so that there are no 
solutions for clusters with a cutoff parameter a > (see Bisnovatyi-Kogan et al. 1998). 
The particular values of a given above characterize the Newtonian curve Mr (a), at T ^ 0. 
In fact, the values aoo{T) and aiim{T), changing for different values of temperature T (see 
center of spirals and maximum values of a for each curve in Fig. 3), in the Newtonian limit 
of T ^ reach the particular values a^o — 2.02 and — 2.87 (see Bisnovatyi-Kogan et 
al. 1998). 

The criteria based on evaluation of entropy and adiabatic invariants are still valid, but 
their application is too comphcated. Analysis made by Bisnovatyi-Kogan et al. (1998) had 
shown that the most convenient criterion of the dynamical stability of rclativistic clusters 
should be based on the investigation of the curves of dependence of the specific binding energy 



Newtonian clusters with another type of cutoff had been studied by Katz (1980). 
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of the cluster Ej, on the central density po, at constant value of Wq. Note that binding energy 
of the relativistic cluster is equal to the total energy of the cluster in Newtonian case, where 
rest mass energy does not appear in the definition of E. 



To analyze the stability of relativistic clusters we need to calculate the specific binding 
energy of the equilibrium models. We have calculated two families of curves, which charac- 
terize dynamical and thermodynamical stability of relativistic clusters with different cutoff 
parameters. The curves -E'{,(po) of specific binding energy at constant temperature T (Fig. 4) 
characterize the thermodynamical stability, while the curves i?f,(po) at constant Wq (Fig. 5) 
give information about dynamical stability of the cluster. Relativistic expression of specific 
binding energy is Ei, = {Nm — M)/Nm, where N is the total number of stars in the cluster 
given by 



and m is the mass of a single star (all stars have the same mass). The number density 
n is expressed as integral in momentum space with similar calculation procedure used for 
obtaining the pressure P and total energy density pc^ (Bisnovatyi-Kogan et al. 1993, 1998). 

The temperature is increasing along each curve in Fig. 5, tending to a finite constant 
value for large values of central density p^. The loss of stability, characterized by the first 
maximum, takes place only for Wq < 15.5. In correspondence of this critical value Wq ^ 15.5, 
the temperature, at large po; reaches a limiting value = 0.0635. This means that no 
dynamic instabilities are present for T < 0.06. At Wq — 16, for example, the limiting 
temperature is equal to Ta — 0.597 and specific binding energy £^b(po) increases monotonously 
until the asymptotic value Eb^a — 0.0312. 

At large po, for models with very large central redshift Zc, there is an asymptotic value 
Eh^a of specific binding energy for each value of Wq. Plotting the function Ei,^a(Wo) from 
Fig. 5 we obtain a more precise boundary of the dynamical stabihty Wo,a = 15.8. Due 
to monotonic dependence of asymptotic (at large po) values of limiting temperature Ta on 
the parameter Wo, similar curve Ei,^a{Ta) from Fig. 4 shows the appearance of dynamically 
unstable clusters at Ta > 0.06. The hmiting curves of specific binding energy Eh^ai^o) and 
Eb,a{Ta) are represented in Figs. 6a and 6b respectively. In the equilibrium configurations 
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with very large central redshift, the temperature is decreasing monotonously with the increase 
of Wo, as may be shown in Fig. 7. Note that the curve T'a(Wo) is approximated with a good 
precision by the power-law relation 

Combining results of numerical investigation of dynamical and thermodynamical stabil- 
ity are represented in Figs. 8, 9 and 10. In Figs. 8a and 8b dynamically stable and unstable 
regions are represented in the planes (T, po) and {T,Zc), respectively. The results plotted 
in the plane {T,Zc) are analogous to ones of the work of Merafina (1999). In Figs. 9a and 
9b thermodynamically stable and unstable regions are represented in the planes (T, po) and 
{T,Zc), respectively. Also in this case the results plotted in the plane {T,Zc) are analogous 
to ones of Merafina (1999). Summary of numerical results on dynamical and thermodynam- 
ical stability analysis is given in Figs. 10a and 10b, where different regions are represented 
in the planes (T, po) and {T,Zc), respectively. It is important to stress the coincidence of 
boundaries between dynamically and thermodynamically stable and unstable configurations 
at large temperatures. Using approximate criteria of dynamical stability we cannot defi- 
nitely judge if these boundaries coincide exactly or there is a small difference between them, 
however the behavior of specific binding energy Ej,/N let us enough confidence in this result. 

Analysis of models with a approaching a^, when maximal densities are expected, had 
shown that they have a very extended halo so, even at very high densities, the cluster 
could remain to be dynamically stable with local Newtonian properties. Similar situation 
should appear when taking into account of degeneracy (Merafina & Ruffini 1990). In the 
Newtonian limit, only nonrelativistic degeneracy is expected, so the situation is not changed 
quahtatively. The logarithmic curve M(po) for a fully degenerate nonrelativistic gas is a 
monotonic line M ~ p^/^, like in a polytropic star with index n — 3/2 and 7 = 5/3, which 
does not show neither thermodynamical nor dynamical instabilities. The thermodynamical 
instabilities on the curve Mr(po) appear only at finite temperatures. 

Consideration of fully degenerate ultrarelativistic particles with equation of state P — qe, 
made by Chavanis (2002b), has sense only for a hmited box because these configurations are 
dynamically unstable with an open outer boundary. The curve £'fe(po) for a given size of the 
box gives results about the onset of dynamical instability in such configurations and has an 
analogy with the corresponding curve at constant Wq. 

In the problem of stellar stability the curve Ms{po), where S is the specific entropy of an 
isentropic star, is used in the static criterion and its maximum determines the boundary of 
dynamical stability of a star (Zel'dovich 1963). The curve M(po) at constant "entropy" was 
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also used for an approximate estimation of the boundary of dynamical stability of truncated 
stellar cluster (Bisnovatyi-Kogan et al. 1993), but such approach is very cumbersome. Using 
the same curve at constant a or is much easier and give very close results. It happens 
however, that at large values of Wq > 11 or at large enough a > 2, the curve M(po) is 
not valid anymore for the estimation of the dynamical stability. The corresponding curve at 
constant a becomes discontinuous and irrelevant. Using the curve Ei,{po) instead of M(po) 
permits to make such estimation in all range of values of the parameter Wq. We have obtained 
numerically that at Wq ^11 maxima of the curves E},{pq) and M{pq) coincide exactly but, 
at larger Wq, maximum of the curve M{po) disappears while the maximum of the curve 
Eh{po) is remaining and shifting to infinite central density at Wo — 15.8. So that all models 
at larger Wq are dynamically stable. Similar relation between Eh{po) and M(po) takes place 
for the curves at constant T. At large values of T > 1 maxima of these curves coincide 
and the curves themselves tend to an asymptotic behavior (see Figs. 2 and 4). At smaller 
temperatures the maximum of Ei){po) is shifting to larger densities (see Fig. 4) and maximum 
of M(po) moves in the opposite side of smaller densities. Following Lynden-Bell & Wood 
(1968) we accept that loss of thermodynamical stability is connected with the maximum of 
the curve Eh^xiPo)- 



6. Conclusions 

The visible similarity between the oscillations of the curve Eb{po) for relativistic clusters 
and the Newtonian isothermal stellar cluster in the box may have a different physical nature: 
in the first case there are two types of oscillations, refiecting the loss of dynamical and 
thermodynamical stability, while the cluster in the box is dynamically stable and oscillations 
are connected only with the onset of thermodynamical instability leading to gravothermal 
catastrophe. 

The analogy between the open cluster with a cutoff in the Maxwellian energy distri- 
bution function and an isothermal cluster in a box is not complete because the first one is 
thermodynamically unstable everywhere while the second one is only after the Antonov's 
point, at density contrast po/Pe ^ 709. Nevertheless, both kinds of oscillations are present 
in the structure of relativistic clusters with a cutoff when we consider clusters with different 
cutoff parameter a and a variable value Wq, which plays the role of the ratio po/Pe (density 
contrast) or — t'l for the cluster in a box. The Figs. 4 and 5 illustrate the behavior of both 
oscillations for the function Eb{po), where oscillations of the curves at constant T add ther- 
modynamically unstable modes and oscillations of the curves at constant Wq give a picture of 
the dynamical stability of the cluster. Calculations show that Newtonian truncated clusters 
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lose their thermodynamical stability at Wq — 6.42, very close to the value —vi — 6.55 for 
clusters in a box. 

We have obtained that only the curve Et,{po) at constant T is valid for the estimation 
of thermodynamical stability, while the corresponding curve M(po) may be misleading. In 
the case of dynamical stability both curves give identical results at lower values of Wq < 11 
while, at larger ones, only the curve £'b(po) may be used. 
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Fig. 1. — Specific binding energy E^/N as a function of Wq, representing the points of 
loss of thermodynamical stabifity of different modes for Newtonian clusters with truncated 
Maxwellian distribution. 



Fig. 2. — Mass M of equilibrium configurations in clusters with a cutoff as a function of 
central density po for different values of temperature T (dotted lines) and cutoff parameter 
a (continuous line). The curves representing dependence M{pq) at constant a have several 
branches at a > 1.5; three branches for a — 2.0 are represented. 




Fig. 3. — Mass M of equilibrium configurations in clusters with a cutoff as a function of 
parameter a for different values of temperature T. 
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Fig. 4. — Specific binding energy Eh/N of equilibrium configurations in clusters with a cutoff 
as a function of central density po ioi different values of temperature T. Each extremum 
corresponds to appearance of new thermodynamically unstable modes. 
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Fig. 5. — Specific binding energy Ei,/N of equilibrium configurations in clusters with a 
cutoff as a function of central density po for different values of parameter Wq. First maxima, 
corresponding to loss of dynamical stability, are present only on curves with Wq < 15.5. 
Each extremum corresponds to appearance of new dynamically unstable modes. 
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Fig. 6. — Specific binding energy Ei,/N oi equilibrium configurations in clusters with a cutoff 
for very large central densities po and central redshifts Zc as a function of Wo (Fig. 6a, left 
side) and (Fig. 6b, right side). The maximum, indicating the loss of dynamical stability, 
corresponds to Wq = 15.8 and T ~ 0.06, respectively. The limiting value of binding energy 
is Eb a = 0.0312. 



-26- 




Fig. 7. — Asymptotic values of temperature T as a function of Wq in the limiting clusters 
with very high central redshifts. 
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Fig. 8. — Regions of dynamical stability and instability in the plane (T, po) (Fig- 8a, left 
side) and in the plane (T, Zc) (Fig. 8b, right side). 
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Fig. 9. — Regions of thermodynamical stability and instability in the plane (T, po) (Fig- 9a, 
left side) and in the plane (T, Zc) (Fig. 9b, right side). 
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Fig. 10. — Regions of dynamical and thermodynamical stability and instability in the plane 
(T, po) (Fig. 10a, left side) and in the plane (T, Zc) analogous to results of Merafina in 1999 
(Fig. 10b, right side). 



